Phase Slips and Interfaces in D— Wave Superconductors. 
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We consider a model (100) interface between two d-wave superconductors. By solving the 
Bogoliubov de Gennes equation on a tight binding lattice, we study the properties of the interface 
as a function of the interface barrier. We contrast the two scenarios: (i) an order parameter phase 
difference of = across the interface, and (ii) a phase slip of S = tt across the interface. We find 
resonant sub-gap structure in the density of states only when there is a phase slip present. We show 
that the local s-wave and "p-wave like" order parameter components are strongly influenced by 
the barrier profile and by the phase slip. The temperature dependence of the local p-wave order 
parameter follows the underlying bulk Tc for p-wave superconductivity implying that this could be 
measured by a suitable tunneling experiment. We also calculate the Josephson critical current as a 
function of the strength of the insulating barrier at the interface. 

Pacs numbers: 74.50. +r, 74.60. Jg, 74.80.-g 



I. INTRODUCTION. 



Over the past few years it has been generally agreed 
that the macroscopic symmetry of the superconducting 
order parameter in high Tc superconductors is unconven- 
tional. The pairing state is almost certainly d-wave [l-p|, 
either pure dx^_y2 or a mixed state such as d+ s |4| or 
d -\- IS P] , which is predominantly d-wave. The d-wave 
pairing state raises many important questions, such as 
the relationship between the pairing state and the physi- 
cal properties of the superconductors, and of course, the 
microscopic origin of the pairing interaction itself. 

In particular it becomes interesting to consider the mi- 
croscopic properties of the interfaces of d-wave symmetry 
systems. This is necessary both to understand the phys- 
ical properties of high Tc materials, and to develop fur- 
ther tests of the pairing state and order parameter. For 
example, it has already been shown that at an normal 
metal to d-wave superconductor [N — D) interface there 
is an extended s-wave component to the order parame- 
ter in the region of the interface and that this survives 
up to a few coherence lengths away from the interface 
P-p|. There is also experimental and theoretical evi- 
dence for microscopic time reversal symmetry breaking 
at such interfaces [|l^,0. The presence or absence of 
such subdominant order parameters can, in principle, be 
used to extract information about the nature of the mi- 
croscopic pairing interaction. For example some models 
will allow both d and s-wave pairing, while other mod- 
els are pair breaking in the s-wave channel. There has 
also been a substantial amount of work concerning the ef- 
fect of d-wave symmetry on the Josephson Effect [[l2|-|l6[ , 
mesoscopic scattering properties oi N — I — D structures 
|l7|-p0| and tunneling density of states for surfaces of d- 
wave superconductors ||2l] , p^ . 

In this paper we examine the self-consistent changes in 
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FIG. 1. A schematic diagram of the system we wish to con- 
sider. Regions I, II and III represent a tight binding lattice, 
where regions II and III are connected to each other via an 
interface (region I). In regions II and III we can define the 
symmetry of the superconducting order parameter to be mis- 
matched by TT. 

the superconducting order parameter near an insulat- 
ing barrier between two d-wave superconductors, i.e. a 
D — I — D interface. The basic geometry is shown in 
Fig. 1. We carry out a detailed study of how the in- 
terface properties evolve as the strength of the barrier is 
increased from a small perturbation to a strong tunnel- 
ing barrier. In this way we can contrast the behaviour 
of interfaces with little or no tunnel barrier, such as twin 
boundaries, to interfaces with a strong tunnel barrier, 
such as tunneling through vacuum or oxygen depleted 
regions at the surface. Our results may also be relevant 
as a model of grain boundary junctions in high Tc mate- 
rials, since Gurevich and Pashitskii have proposed that 
these junctions are essentially D — I — D junctions whose 
impurity barrier is simply dependent up the angle of the 
grain boundary ||23| . 

In our calculations we will concern ourselves with two 



scenarios (i) where there is no phase difference in the or- 
der parameter over the interface {6 = 0) and (ii) where 
there is an imposed phase shp of vr over the interface 
{9 = vr), as illustrated in Fig. |]. For both of these sce- 
narios we focus our attention on the magnitude of the 
extended s-wave and "p-wave like" components of the 
order parameter in the region of the interface. There is 
already a great interest in the role extended s-wave com- 
ponents have to play in the region of twin-boundaries, for 
example W. Belzig, C. Bruder and M. Sigrist [O have 
studied how the local density of states is modified, as one 
moves across a twin boundary when the bulk order pa- 
rameter is d ±is and the effect of temperature on these 
calculations. 

The plan of this paper is as follows. In Sec. II we de- 
scribe briefly the model Hamiltonian, and the methodol- 
ogy used to solve the Bogoliubov de Gennes equation in 
our system. Then in Sec. Ill we will describe, in detail, 
the system of interest and define the different compo- 
nents of the superconducting order parameter which will 
be relevant. In Sec. IV we present results for a perfect 
interface, i.e. the limit that the tunnel barrier height is 
zero. We note that in this clean case an order parame- 
ter phase change oi 9 = tt across the "interface" corre- 
sponds to a phase slip of tt in a perfect superconductor. 
We establish that this phase slip is a stable soliton like 
solution of the Bogoliubov de Gennes equation and we 
obtain the local density of electronic states in the vicin- 
ity of the phase slip. We then proceed to analyze the 
different components to the superconducting order pa- 
rameter at the centre of the phase slip and how these 
change as we change the temperature of the system. In 
Sec V we turn our attention to D — I — D interfaces, 
calculating the properties of the interface as a function 
of the tunnel barrier height for both the scenarios 9 = 
and 9 = TT. Finally we conclude with comparisons of the 
local particle density of states in the region of the inter- 
face. Using these local particle density of states we also 
estimate how the Josephson critical current is changed 
as the tunnel barrier strength is changed. 



II. THE BOGOLIUBOV DE GENNES EQUATION. 

Our starting point is the Bogohubov de Gennes (BdG) 
equation. 



where Hij is the normal part of the tight binding Hamil- 
tonian and Aij is the anomalous part, given by the BCS 
gap function. 

In the case of the cuprates we have dx2-y2 pairing on 
a square two-dimensional lattice. The simplest model 
Hamiltonian which leads to this pairing state is the non- 
local attractive Hubbard model, 



^ = E' 
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where, as usual, c\^ and Cia are the electron creation and 
annihilation operators at site i for spin a and fii is the 
number operator at site i. For simplicity we shall assume 
that both tij and Uij are limited to nearest neighbors 
only. We shall neglect any on-site interactions Uu. Also, 
for simplicity, we neglect retardation in the attractive 
interaction. 

The BdG equation, Eq. |l| is easily derived as the ap- 
propriate mean field factorization of the Hubbard model 
Eq. H The normal part of the Hamiltonian Eq. |^, H^ , 
is given by 
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where /i is the chemical potential. The on-site energy, e^, 
includes the bare on-site energy and the normal Hartree 
potential 



/ , ^ij'^^jj 



(5) 



where rijj is the charge density at site j. Similarly the 
self-consistent hopping t^ ■ is the bare hopping plus a 
Hatree-Fock exchange term 



'^ij ~ ''ij + n'-^ ij''^ij ■ 



(6) 



The charge densities nn and n^ are determined from the 
eigenvectors of the Hamiltonian Eq. n^ and are given by 
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where u" and w" are the particle and hole amplitudes on 
site i associated with eigenenergy En- The Hamiltonian 
is 



-"y - I A*- -H* 
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(2) 



where the sum is over all of the negative eigenvalues En 
only. 

The off-diagonal part of the mean field Hamiltonian 
Eq. is the pairing potential or gap function, A^ . Self- 
consistently this is determined from 



Ay — Uij Tij 



(8) 



where the anomalous density is 
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and again the sum only includes eigenvalues En up to the 
condensate chemical potential (/x). 

In the case of a bulk square lattice with nearest neigh- 
bour attractions C/y and nearest neighbour hopping the 
solutions to the above set of self-consistent equations are 
well known. Depending on the band filling and temper- 
ature the gap function A^ may be pure d^^-y^, pure 
extended s-wave ot a d + is mixed state [Q . In this pa- 
per we shall work exclusively near half filling where only 
the pure d-wave state is stable. 

To solve the above system of self-consistent equations 
we employ the Recursion Method ||^,|2^,^,g^ following 
exactly the method outlined in |^ . The advantage of this 
method is that it is a purely real-space calculation, and 
so it is not necessary to assume periodic boundary condi- 
tions. For example, in the interface geometry studied in 
this paper we have two semi-infinite lattices joined at an 
interface. The recursion method also gives directly the 
local density of states, N{E), at any given site, which, as 
we shall show below, assists in the physical interpretation 
of the self-consistent numerical results. 



III. INTERFACE GEOMETRY. 

In the cuprate superconductors there are a number 
of important intrinsic weak links, including twin bound- 
aries, low angle grain boundaries, and edge dislocations. 
There are also a number of man made interfaces designed 
for specific applications, including break junctions, high 
angle grain boundary junctions, ramp junctions, and 
junctions with electron and ion beam irradiated barrier 
regions. In separate papers we have investigate simple 
S — N, D — N and S — D interfaces [^ , and examined in 
detail the behaviour of a I? — 13 high angle grain bound- 
ary junction [ P8[ . In this paper we seek to examine how 
the interface properties are modified by the presence of 
an insulating barrier. For example the insulating bar- 
rier may be caused by oxygen depletion in the region of 
the interface, as suggested by Gurevich and Pashitskii in 
their model of grain boundaries [E3[ . 

With this motivation, we shall study in detail the sys- 
tem shown in Fig. 0(a). We have two semi- infinite d- 
wave superconducting regions (II) and (III) which arc 
connected to each other via an interface (I). We study 
(100) oriented interfaces as shown in Fig. P(a). In the in- 
terface region (I) we define the strength of the barrier, V , 
by setting e[ — V 'm this region to be finite, compared to 
the rest of the system where e'j = 0. This barrier height 
could physically represent the amount of oxygen deple- 
tion at the interface, or the vacuum or insulating barrier 
width in the case of (100) surface tunnel junctions. We 



shall examine the changes in the interface properties as 
the barrier height V is varied from zero to a large value. 

We shall study the Josephson energy of the junction by 
comparing the cases of either zero or tt phase difference 
between the bulk order parameters far from the inter- 
face. Fig. |l|(b) illustrates the case of a tt phase difference 
in the d-wave order parameter. We carry out the cal- 
culation by imposing a fixed bulk d-wave order parame- 
ter Ajj far from the interface, and then self-consistently 
computing the A^ in the region near to the interface. 
Typically we perform the self-consistent calculations in 
a region of width five coherence lengths on either side of 
the junction, confirming that after this distance the self- 
consistent parameters have returned to their bulk values. 

The self-consistent order parameters Ay are computed 
for each nearest neighbour bond < ij > in the lattice 
(Ay — Aji). However interpreting these results is sim- 
plified if we decompose these Ay into components of dif- 
ferent symmetries. For this reason we define the following 
four linear combinations at each site i, 
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Here the sites ji . . . J4 are the four nearest neighbours of 
site i, say counting clockwise. The first two combinations 
are simply dx2-y2 and extended s-wave pairing respec- 
tively. The second two correspond to components of the 
Ay which are odd about site i. For this reason they may 
be said to be "p-wave like" components of the order pa- 
rameter. It is important to note that these components 
do not truly correspond to p-wave pairing since they are 
spin singlet not triplet and there is no spontaneous break- 
ing of spin rotational symmetry. These two terms are 
more closely related to gradients of s or d-wave order pa- 
rameters. They also have some relationship to Yang's rj 
pairing order parameter p9| which is a spin singlet order 
parameter, but one that alternates sign on neighboring 
lattice sites. This has also been referred to elsewhere as 
antif err o superconductivity |pQ,pl|. Here A'^^^ and A^^"' 
correspond to an alternating non-local order parameter 
Ay, unlike //-pairing which corresponds to a staggered 
Aii on-site gap function. Note that because the system 
in Fig. ^a) is perfectly periodic in y we shall always have 

^!f''^ = 0. 

For our chosen model parameters only the d-wave com- 
ponent A^ is non-zero in the bulk far from the inter- 
face. However components of the other types can arise 
near the interface driven by the local symmetry break- 
ing. The self-consistent order parameter can also have 



spontaneous symmetry breaking phase transitions as a 
function of temperature or chemical potential, for exam- 
ple breaking of time reversal symmetry into an d+is type 
state. 



IV. INTERFACE PROPERTIES: NO BARRIER 



Having outlined in previous sections the method for 
calculating and analyzing the properties of the order 
parameter, we now consider the most simple situation: 
V = throughout the system. In this case there is no 
barrier at all at the interface /. If the bulk phase differ- 
ence across the system 9 is zero then we obviously have 
just a simple bulk system everywhere with the bulk d- 
wave order parameter. 

On the other hand if 6* = tt (Fig. |l](b)) then we have a 
phase slip in the unit cell. This is the superconduct- 
ing analogue of a Bloch wall separating magnetic do- 
mains. By solving the Bogoliubov de Gennes equations 
self-consistently we can calculate the width of the phase 
slip region and the variation of the order parameter across 
the phase slip. 

Fig. ^ shows our self-consistently determined order pa- 
rameters in this phase slip geometry. Since all quanti- 
ties only vary in the x direction, we can plot |A''''(2;)|, 
A(^)(x)| \A^p^^{x)\ in the region of the interface {x = 
0), noting that \A^Py\x)\ = 0. We started the self- 
consistent calculation by imposing a d-wave order pa- 
rameter throughout the system, but with a phase change 
of 77 at a; = 0. The solution was then determined self- 
consistently by iterating the BdG equations starting from 
this state. We find that the phase slip is a stable self- 
consistent solution of the BdG equations. 

In Fig. g we can see how the | A^'') {x) \ (solid line) varies 
in the region of the interface. It is clearly a minimum at 
the interface and quickly reaches it's bulk value, oscillat- 



ing as it does |3^,Q. The extended s-wave component 
A^*''(x)| (insert) is also a minimum at the interface, but 
it then reaches a maximum close to the interface and then 
reduces to zero in the bulk. On the other hand |A(P=")(a;) 
(dashed line) has a strong maximum at the interface. 




FIG. 2. Profiles of different symmetry components of 
the superconducting order parameter in the presence of 
a e = -K phase slip, |A('^'(a;)|/|A('*)(oo)| (solid line) 
and |A'*'=">(x)|/|A''''(oo)| (dashed line). The insert shows 
I A*"" (x) I in the region of the interface. For all of these figures 
T = and the barrier height is V = 0. 




FIG. 3. The effective hopping t'j as a function of distance 
from the interface {x = 0). The hoppings are t'ji{x) (solid 
line), t'j2{x) (dotted line), t'j^{x) (dashed line) and t'j4^{x) (dot- 
ted fine), where j'l to J4 correspond to the four different near- 
est neighbour bonds. All parameters are the same as in Fig. 
2. 

It is also interesting to consider the self-consistent hop- 
ping, ij -, in the region of the barrier. Again, considering 
that each site i has four neighbours ji . . . J4 we can plot 
the four components t'^^^ . . . t'^j^ as a function of position. 
Changing notation slightly we can think of these as four 
functions of position, x, t'Ax) {j — I. .4). In Fig. 3 the 
dashed and solid lines show t[{x) and t'^{x) respectively 
(the hopping perpendicular to the interface) and the dot- 
ted hue is t2{x) {— t'^{x)), (the hopping parallel to the 
interface) . The plot shows that in the region of the in- 
terface the hopping is modified from its bulk value. It 
is this modification of the hopping which helps stabilize 
the phase slip, this is analogous to the case of conven- 
tional local s-wave pairing where modification of the lo- 
cal density in the region of the phase slip enhances its 
stability [p3. The oscillations in the ij shown in Fig. 3 
simply correspond to Friedel oscillations in the non-local 
Hartree-Fock exchange term ( '^J*'^ ) . These oscillations 
are exactly analogous to the Friedel oscillations around 
an impurity or surface in a normal system, even though 
there is no normal barrier V here. 

The results shown in Figs. and M were performed 
with T <^ T^, the bulk d-wave transition temperature. 
We now proceed to examine the temperature depen- 
dence. In Fig. 4 we have plotted the normalized temper- 
ature dependence of |A(p-)(0)|2 (circles) and |A('*)(oo)|2 
(crosses), these being the p-wave like order parameter at 
the interface and the bulk d-wave order parameter re- 
spectively. As we can see |A(p^)(0)P is a linear function 
of temperature over a wide range (from T ^ O.lt to 0.4t). 
Extrapolating this linear dependence we can state that 
|A(p=)(0)P would go to zero at a temperature of around 
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FIG. 4. The temperature dependence of the p-wave com- 
ponent of the order parameter, ' ,p_^. i — (circles), and the 



d-wave component, 
units of t. 
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(crosses), T is measured in 



TP « 0.5t, if it were not coupled to the d-wave order 
parameter. From Fig. 4 we infer from the fact that 
|A(p=)(0)P deviates from a straight hne above T k, 0.5i 
that it is enhanced above T-P by couphng to the d-wave 
order parameter. 

This unusual temperature dependence of the p-wave 
hke component can be explained by examining the phase 
diagram for a perfect square tight binding lattice. It is 
known that at half filling we have a ground state with a 
stable d-wave order parameter. However one can also de- 
fine a separate Tc for extended s-wave or p-wave states. 
Near half filling T^ is greater than T^ or T^ and so only 
the d-wave state occurs in the bulk [|4J . However if the 
d-wave and extended s-wave order parameters were sup- 
pressed somehow then the system would in principle go 
superconducting at T^. Interestingly the T^ for the "p- 
wave like" singlet order parameter A^^-* is the same as 
that of true triplet p-wave order parameter, which fol- 
lows from the identical form of the gap equation in both 
cases [|4|. In effect our phase slip has forced A'^''-' and 
A*^") to vanish at a; = 0. This allows the intrinsic un- 
derlying p-wave state to become apparent. It would be 
interesting to test this underlying p-wave state experi- 
mentally, for example by tunneling experiments to see 
whether or not there is a finite gap inside a supercon- 
ducting phase slip. If there is such a gap it would show 
that the pairing mechanism has attractive components 
for p-wave pairing as well as d-wave. This would only 
be true for some specific model pairing interactions, such 
as our nearest neighbour Hubbard model. 



V. INTERFACE PROPERTIES: DEPENDENCE 
ON BARRIER. 

Having considered the most simple case (e^ = for all 
sites) we now focus our attention on what happens if we 



have a finite tunnel barrier {V ^ 0) in region (I) of our 
system (Fig. |l|). In this case it is also interesting to con- 
sider the situation where there is no phase slip over the 
interface and compare this to the result obtained when 
there is a phase slip. In this section we study three dif- 
ferent strengths of barrier, where in region (I) we have 
defined V/t= 1,2,5. 
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FIG. 5. Profiles of different symmetry components of the 
superconducting order parameter in the presence of no phase 
slip, ^ = 0, (solid lines), and a 9 = n phase slip (dashed 
lines). The upper graph plots |A*P"'(a;)|/|A''*^(oo)| and the 
lower graph plots |A* '(a;)|/|A'''^(cx))|. The barrier interface 
height is V = t and temperature is T = 0. 

In Fig. g we compare both |A(p="'(x)| and |A('')(a;) 
for the cases of (i) there is a phase slip across the in- 
terface (solid line) and (ii) when there is no phase slip 
(dashed line), for V/t = 1. From this figure we can see 
that as in the case of no barrier the phase slip means 
that |A('')(0)| = 0. This then, effectively, induces a fi- 
nite |A(p="^(0)| at the interface. However when there is 
no phase slip (dashed lines) we can see that although the 
A('')(a;)| is diminished in the region of the interface it 
does not go to zero, conversely now, due to the symme- 
try of the problem, |A(p=='(0)| must be zero. We have 
also calculated the extended s-wave component to the 
order parameter in the region of the interface. We find 
that the extended s-wave component is a maximum, at 
the interface, when there is no phase slip and zero when 
there is a phase slip, similar to Fig. 0. 

Figs. ^ and show the different components to the or- 
der parameter for larger strengths of barrier {V/t = 2,5). 
What we see in these two figures is that as the barrier 
height is increased the differences between the 9 — 
and — t: cases becomes smaller. In the limit of a very 
high strength impurity barrier. Fig. u[ there is essentially 
no difference between the order parameter components 
whether a phase slip is present or not. Fig. O, however, 
is more interesting. We see that when there is no phase 
slip |A^'^)(a;)| is reduced in the region of the barrier, and 
as expected |A^p^)(0)| is zero. In the presence of a phase 




FIG. 6. Profiles of different symmetry components of tlie 
superconducting order parameter, for an intermediate barrier 
lieiglit, V = 2t. The plotted curves are otherwise exactly as 
in Fig. 5. 




FIG. 7. Profiles of different symmetry components of the 
superconducting order parameter, for a large barrier height, 
V = 5i. The plotted curves are otherwise exactly as in Fig. 
5. 

slip, as in all the other cases, the jA^'^-'fO)! is zero at 
the phase slip. But, unlike Figs. § and ||, |A(p^)(0)| is 
also zero. This shows that we have gone through three 
regimes at the interface. These being (i), low scatter- 
ing at the interface where, |A('')(0)| is finite when there 
is no phase slip and |A'p^)(0)| is finite in the presence 
of a phase slip, (ii), intermediate scattering at the inter- 
face where, |A'^''^(0)| is still finite when there is no phase 
slip but now |A(p^)(0)| is zero in the presence of a phase 
slip and (iii), large scattering at the interface, |A('')(0)| 
is now zero when there is no phase slip and |A'^^="^(0)| is 
also zero in the presence of a phase slip. 

To see the physical origin of these three regimes we 
study the local particle density of states, N{E), in the re- 
gion of the interface. Figs. 0, M, 00 and nTlshow the den- 
sity of states N{E) for four strengths of barrier, V/t — 0, 



1, 2 and 5, respectively. In these four figures we have 
plotted both the density of states in the presence of a 
phase slip, N'^{E), (thick solid line) and in the absence 
of a phase slip, N^{E), (thin solid line). We also show 
the densities of states for four values of x, the distance 
from the interface. 
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FIG. 8. Local particle density of states in the region of the 
interface when the barrier height is V = 0. We consider both 
a phase slip (thick solid line) and no phase slip (thin solid 
line). The curves plotted are (a) a; = 0, (b) a; = 1, (c) a; = 2 
and (d) x — 3. 



In the first case we consider V/t = (Fig. g) . The first 
thing to note is that the density of states in the absence of 
a phase slip corresponds to the pure bulk lattice d-wave 
density of states. This is simply because ioi V = and 
6 — there is no interface. However the density of states 
in the presence of a phase slip, N^{E) does change in as 
we move away from the interface. Comparing N^{E) and 
N'^{E) we see a dramatic shift in the density of states in 
the region of the Fermi energy. We can see in Fig. S(a) 
that there is a strong shift of spectral weight from the 
peaks into the region near the Fermi level. We see that 
the gradient of the local particle density of states, in the 
region of i? = is dramatically changed, and the size of 
the gap, (the distance between the two peaks) is strongly 
reduced in the presence of a phase slip. Moving away 
from a; = (Figs 8(b), 8(c) and 8(d)) we see that these 
large changes in density of states gradually diminish with 
X until N^{E) and N'^{E) become nearly identical. In- 
terestingly the energy region near .E = is most strongly 
affected when x is large, corresponding to the long-range 
nature of the energy states associated with the d-wave 
gap node [p5 36|. 



There is a very strong contrast between the case T^ = 
in Fig. g and the case V/t = 1 shown in Fig. pi In Fig. O 
we see that N^{E) (thin solid line) has no subgap struc- 
ture at all, and this is quite different to the dramatic 
subgap structure obtained for N'^{E). 
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FIG. 9. Local particle density of states in the region of the 
interface when the barrier height is weak, V = t. The plotted 
curves are otherwise exactly as in Fig. 8. 
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FIG. 10. Local particle density of states in the region of 
the interface when the barrier height is intermediate, V — 2t. 
We consider both a phase slip (thick solid line) and no phase 
slip (thin solid line). The curves plotted are (a) x = 1, (b) 
a; = 2, (c) a: = 3 and (d) x = 4. 

As we increase the barrier strength to V/t = 2 ( 
Fig. |o|) we still see that there is no subgap structure in 
N^{E) whereas N'^{E) still has some subgap structure, 
although this structure is not as dramatic as in Fig. ^. 
In Fig. |l^ we have only plotted N°{E) and N^'iE) for 
x 7^ because the x — density of states is dominated 
by states well away from the gap region. 

In Fig. ^ we see that in the presence of a very strong 
barrier, V/t = 5, N°{E) = N^'iE) suggesting that the 
two sides to the interface have become effectively decou- 
pled. Also in this limit there is no resonant structure in 
either N"{E) or N'^iE). 

Finally we estimate the Josephson critical current as a 
function of the barrier height, by determining the Joseph- 
son energy. This is defined by 



LU 




oj N^V^-^ 



E/t 

FIG. 11. Local particle density of states in the region of 
the interface when the barrier height is strong, V = 5t. The 
plotted curves are otherwise exactly as in Fig. 10. 



Ej = l {E{7r) - E{0)) 



(10) 



where we estimate the energy difference from the inte- 
grated densities of states 



E{<t>) 



N{(j)){E)dE 



(11) 



for (f) — n and cj) — respectively. The critical current is 
approximately proportional to this energy difference. 



Ic oc Ej 



(12) 



where the constant of proportionality depends on the 
form of the £'(</>) relationship. In the simple Josephson 
tunneling case E{(j)) — E{0) ~ Ej sin (j), and Ic — Ej, but 
in general other forms of E{(l)) are possible. For example 
we have found that E{(j)) is approximately sawtooth in 
grain boundary junctions ||28|1 and there is experimental 
evidence for non-sine wave behaviour ||37|1 . 

In Fig. |lj we have plotted the calculated critical cur- 
rent (Ic) as a function of the strength of the impurity 
potential (e/t) in region /. In this figure we have nor- 
malized the critical current with respect to the critical 
current when no barrier is present. What we again ob- 
serve is consistent with the three different regimes (i) , lit- 
tle or no normal scattering at the interface, where there 
is no sub gab resonant structure but N^{E) ^ N^{E) 
for |_E| < lA^"^)], (ii), moderate normal scattering at the 
interface, where there is sub gap resonant structure in 
N'^{E) but not in N^{E) and finally (iii), large normal 
scattering at the interface, where N^{E) « N'^{E). From 
this it is possible to deduce that, as one would expect for 
a high strength interface barriers, the Josephson Critical 
Current approaches zero. The three regimes are visible 
in the calculated Ic of Fig. ^3: in the first weak decrease 
in Ic with barrier height, the intermediate regime where 
there is a strong decrease in Ic, and in the strong scat- 
tering regime where Ic —> 0- 
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FIG. 12. The normalized Josephson critical current as a 
function of interface barrier height, V/t. 



VI. CONCLUSIONS. 

In this paper we have investigated the properties of 
a fi-wave order parameter in the region of a (100) in- 
terface between two d^2_y2 superconductors. We have 
determined how the barrier height at the interface influ- 
ences the local density of states and Josephson current. 
In particular we could identify three distinct regimes with 
qualitatively different behaviour, corresponding to weak, 
intermediate and strong barriers, respectively. 

We have also calculated the subdominant extended s- 
wave, and p-wave like order parameter components at the 
interface. In particular we have shown that the p-wave 
component has a temperature dependence which reveals 
the underlying bulk Tc for p-wave pairing, even though 
this is normally hidden by the higher bulk d-wave Tc- 
In principle a tunneling experiment could be carried out 
to test whether or not this p~wave component occurs in 
real systems. Such an experiment would provide qualita- 
tively new information about the form of the microscopic 
pairing interaction. 
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